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ABSTRACT 


In  this  paper  we^considerijlinear  multiple  objective  programs  with 
coefficients  of  the  criteria  given  by  intervals.  This  class  of  problems 
is  of  practical  interest  since  in  many  instances  it  is  difficult  to  deter- 
mine precisely  the  coefficients  of  the  objective  functions.  A subproblem 
to  test  if  a feasible  extreme  point  is  efficient  in  the  problem  considered 
is  obtained.  A branch  and  bound  algorithm  to  solve  the  subproblem  as 
well  as  computational  results  are  provided.  Extensions  are  discussed. 


LINEAR  MULTIPLE  OBJECTIVE  PROBLEMS  WITH  INTERVAL  COEFFICIENTS 
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1.  Introduction 

The  linear  multiple  objective  problem  (LMOP)  consists  In  choosing  alter- 
natives from  a polyhedral  set  considering  simultaneously  conflicting  linear 
criteria.  The  LMOP  Is  vnrltten  as 

(LMOP):  Max{Cx;xEF} 

where  F » {xeR™  : Ax*b,  x^O) , C and  A are  p by  n and  m by  n matrices  respec- 
tively, and  beR®.  A point  x'eF  is  said  to  be  efficient  in  (LMOP),  or  equiva- 
lently with  respect  to  C,  if  there  is  no  xeF  such  that  Cx  > Cx“  with  at  least 
one  strict  inequality.  The  set  of  efficient  points  is  denoted  by  EF(C)  and  is 
considered,  in  this  paper,  to  be  the  solution  set  to  (LMOP). 

The  criteria  as  given  by  C can  be  seen  either  as  the  linear  utilities  of 
p decision  makers  or  as  the  p object.ve  functions  of  a single  decision  maker. 
Several  real  and  potential  applications  of  multiple  criteria  problems  have  been 
recently  reported  In  areas  such  as  water  resource  planning  and  facility  location 
[8],  scheduling  of  nursing  resources  [9],  zero  defects  program  [2],  evaluation 
of  urban  policy  [1],  investment  decision  making  [18],  resources  allocation 
[21],  energy  planning  [38],  macroeconomic  policy  [34],  forest  management 
[33],  location  of  public  facilities  [28],  activity  planning  [26],  and 
corporate  financial  management  [24]. 

Different  approaches  have  been  suggested  to  solve  decision  problems  when 
more  than  one  objective  is  considered.  Survey  articles  by  Roy  [29],  MacCrlmmon 
[25],  Kombluth  [23],  Gal  [15],  Starr  and  Zeleny  [31],  Flshbum  [14],  Farquhar 
[13]  and  recent  books  by  Keeney  and  Raiffa  [22],  Cohon  [8],  Cochrane  and 
Zeleny  [7],  Zlonts  [37],  Starr  and  Zeleny  [30],  and  Hwang  and  Masud  [19] 
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present  an  extensive  coverage  of  the  methods  proposed.  In  this  paper  we  con- 
centrate on  the  approach  that  defines  the  set  of  efficient  points  as  the  solu- 
tion set  of  the  multiple  criteria  problem.  In  particular,  we  are  Interested 
in  the  set  of  efficient  extreme  points  of  F that  we  denote  by  EF^^(C). 

Algorithms  devised  by  Charnes  and  Cooper  [6],  Ecker  and  Kouada  [11],  Evans 
and  Steuer  [12],  Isermann  [20],  Gal  [16],  and  Yu  and  Zeleny  [36]  can  be  used 

to  obtain  EF  (C) . 

ex 

In  practical  applications  it  is  usually  difficult  to  determine  the  coeffi- 
cients of  the  criteria  matrix  C because  they  are  either  specified  subjectively 
by  the  decision  maker (s)  or  they  are  obtained  through  procedures  requiring 
subjective  answers  to  questions  posed  by  the  analyst.  In  other  Instances, 
the  criteria  are  obtained  by  least  square  minimization  or  by  linear  regression. 

An  important  example  where  imprecision  of  coefficients  is  known  to  exist  and 
where  the  consideration  of  multiple  criteria  is  vital,  is  public  sector 
decision  making.  It  Is  very  difficult  to  estimate  premiums  applied 
to  market  prices  to  obtain  social  costs,  shadow  price  of  investment  and  rates 
of  saving  of  different  segments  of  the  society  ([5]  and  [10]).  The  effect 
of  the  inaccuracy  in  the  elements  of  C is  usually  hard  to  evaluate  and  no 
computationally  effective  method  is  available  to  perform  a full  parametric 
analysis  in  (LMOP) . 

In  this  paper  we  propose  the  use  of  Interval  estimates  for  the  elements  ® 
of  C instead  of  the  current  practice  of  point  estimates.  In  most  Instances 
decision  makers  may  feel  more  comfortable  in  specifying  intervals  than 
points;  furthermore,  information  available  from  a statistical  analysis 
for  determining  the  criteria  may  be  better  utilized  by  giving  the  coefficients 
in  the  form  of  Intervals.  These  can  be  seen  as  confidence  intervals 
for  the  c^j's.  In  fact,  situations  where  the  proposed  approach  is  of  parti- 
cular Interest  is  when  the  elements  of  0 are  known  to  be  stochastic. 
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The  discussion  above  suggests  the  consideration  of  the  following  linear 
multiple  objective  problem  with  interval  coefficients: 


(P): 


Max{Cx  : xeF,  Ce'I'} 


where  F is  as  previously  defined  and  is  the  set  of  p by  n matrices  with 

elements  c^^^  in  the  interval  f j » j 1 and  j”l,2,...,n.  The 

lower  and  upper  bounds  and  are  given  real  numbers.  Problem  (P)  represents 

a family  of  (LMOP)'s,  with  one  problem  for  each  C in  The  solution  set  to 

(P)  is  defined  as  the  set  of  points  in  F efficient  with  respect  to  every  Ce4’ 

and  is  denoted  by  EF.  VJe  also  refer  to  it  as  the  set  of  efficient  points  in 

(P) . The  definition  of  EF  Implies  that  EF  = H EF(C). 

Ce$ 

The  plan  of  the  paper  is  as  follows.  In  section  2 properties  of  (P), 
Including  connectedness  of  EF  and  the  existence  of  efficient  points,  a'te 
discussed.  In  section  3 the  case  where  all  lower  bounds  are  nonnegative 
is  considered.  An  algorithm  to  determine  EF^^  is  obtained  in  section  4.  It 
is  shown  that  to  determine  the  set  of  efficient  extreme  points  in  (P)  it  suffices 
to  use  any  of  the  algorithms  to  solve  (LMOP) , for  a fixed  C,  in  conjunction  with 
a subroutine  that  solves  a test  problem  at  each  extreme  point.  Although  the 
test  problem  is  nonlinear,  an  effective  implicit  enumeration  scheme  that  only 
considers  linear  programs  has  been  developed  to  solve  it.  Computational 
results,  extensions,  and  topics  for  further  research  are  discussed  in  section  5. 
Straightforward  proofs  are  omitted  for  the  benefit  of  compactness  of  the  paper. 

The  following  notation  is  given  for  future  references.  Lower  case 
letters  are  used  to  denote  vectors.  Superscripts  differentiate  vectors 
and  subscripts  indicate  the  components  of  a vector.  The  partial  ordering 
relation  x ^ y means  Xj  i with  at  least  one  strict 
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inequallty.  Matrices  are  denoted  by  capital  letters.  The  row  and 
column  of  a matrix  C are  written  as  and  respectively.  Given 

two  matrices  C*  = C*  “ notation  C’  ^ C ^ C*  indicates  that 

the  matrix  C “ has  each  element  c^^  satisfying  cjj  ^ c^j  < c* ^ . 

Given  a set  S,  we  denote  the  subset  of  its  extreme  points  by 
Max{....}  (max{...})  represents  a multiple  (a  single)  objective  maximization 
problem.  Consider  the  problem  Max{Cx  : Ax=b,  x^O}  and  let  x®  be  an  extreme 
point  of  the  constraint  set.  Without  loss  of  generality  the  basis  B, 
corresponding  to  x®,  is  assumed  to  be  composed  by  the  first  columns  of  A,  x® 
is  written  as  x®  = (x®  ,x®  ),  A and  C are  accordingly  partitioned  as 

A ” (B,A^)  and  C = (Cg.Cjj)  and  the  linear  system  multiplied  by  B ^ becomes 

T oB  . tN  .N  r . tN  „-1.N  , r -1. 

I X + A X = b where  A = B A and  b = B b. 


2.  Properties  of  (P) 

The  set  $ is  contained  in  the  vector  space  of  p by  n matrices  and  is 
characterized  as  follows. 


Proposition  2.1;  <li  is  a convex  set  with  extreme  points  having  each  element 
Cj^j  at  its  upper  or  lower  bound. 

For  each  Ceit  let  K(C)  = {pcR”  : Cp^O}  and  let  K('t')  *=  U K(C).  k(C)  and 


Ce$ 

K($)  are  called  the  preference  cones  associated  with  matrix  C and  problem  (P) , 
respectively.  Note  that  the  elements  of  K(C)  are  the  directions  of  preference, 
l.e.,  if  y “ p+x  for  some  peP(C)  then  Cy  ^ Cx.  Preference  cones  play  an 
important  role  in  general  multiple  criteria  optimization  (see  [35]  and  [4]). 

They  provide  an  Insightful  geometrical  Interpretation  and  are  intimately  related 
to  the  connectedness  and  existence  of  efficient  points.  Denote  by  M the  subset 
of  '"atrices  of  t having  all  elements  of  each  column  at  the  upper  bound  or  at 
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the  lower  bound.  Hence,  if  CeM,  for  j=l,2,...,n  either  or 

where  U [y^jl  and- L = The  maximum  number  of  elements  in  M is  2". 

Proposition  2.2;  KC*?)  = U K(C)  A K(M)  . 

CeM 

Proof:  Assume  peK($).  Thus,  there  is  Ce^  such  that  Cp  0.  Consider  the 

matrix  c’  = defined  for  j=l,2 n as  c'j  = if  p^  < 0 and  = U^j 

otherwise.  Then,  Cp  ^ Cp  ^ 0 and  peK(c*).  Since  by  definition  c'eM  it  follows 
that  K(il>)  C K(M)  . The  reverse  Inclusion  relation  is  a direct  consequence  of 
the  definitions  of  K(<J))  and  K(M). 


Proposition  2.2  shows  that  from  the  infinite  family  of  problems  represented 
by  (P)  we  just  have  to  consider  the  finite  family  of  (LMOP)'s  having  C in  M. 


Corollary  2.3:  EF  - H EF(C)  A EF(M) . 

CeM 

Proof:  By  definition  EF  « O EF(C).  Since  MCt",  EFCEF(M).  Assume  xeF  and 

Celt  “ 

x^EF.  Thus,  [{x}  + K(iIi)]nF  ^ {x}.  By  Proposition  2.2  it  follows  that 
[{x}  + K(M) ] OF  y {x}  and  therefore  x^EF(C)  for  some  CeM.  Consequently 
EF(M)CEF  and  hence  EF  = EF(M). 

Proposition  2.2  and  Corollary  2.3  imply  that  determining  the  efficient 
points  in  (P)  is  equivalent  to  solve 

(P'):  Max{Cx  : xeF,  CeM}. 

This  last  problem  represents  a finite  but  large  family  of  (LMOP)’s.  The 
algorithm  developed  in  section  4 solves  (P')  using  an  efficient  implicit 
enumeration  scheme. 

The  linear  multiple  objective  ,jroblem,  LMOP,  has  been  extensively  studied, 
(e.g.,  [12] , [27] , [16] , [36] , [ 39] , and  [11]).  It  is  well  known  that  EF(C) 

is  closed  and  connected,  the  existence  of  an  efficient  point  implies  the 
existence  of  an  efficient  extreme  point,  the  preference  cone  K(C)  is  convex 


X. 
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and  whenever  a point  In  the  relative  Interior  of  a face  Is  efficient  the  whole 
face  also  Is.  The  following  proposition  establishes  the  corresponding  results 
in  problem  (P) . 


Proposition  2.4;  a)  EF  is  closed, 

b)  If  a point  in  the  relative  interior  of  a face  of  F is 
efficient,  the  whole  face  also  Is, 

c)  If  EF  ^ <J),  there  is  an  extreme  point  of  F efficient  in 


d)  is  not  necessarily  convex  and, 

e)  EF  Is  not  necessarily  connected. 

Proof:  a),  b) , and  c)  follow  trivially  from  the  properties  of  (LMOP)  and 
Corollary  2.3.  d)  and  e)  are  illustrated  with  examples.  Let  4>  be  In  the 
space  of  1 by  2 matrices  and  let  = 1,  = .5,  *^22  “ 

Consider  p*  = (.6,-1),  C’  = [1, .5]e$,  p*  = (-1.8,1),  and  C^  = t-S.lle^. 

Then,  C'p*  > 0,  C^p^  > 0 but  p •=  p*+p*  = (-1. 2,0)^K(iJ>) . It  Is  not  difficult 
to  construct  examples  where  EF  is  not  connected  when  the  closure  of  K(4>)  contains 
a non-trlvlal  si’bspace  as  is  shown,  for  general  multiple  objective  optimlia- 
tlon  problems,  in  [4].  However  in  (P) , even  when  this  is  not  the  case,  EF  may 


be  not  connected.  For  example,  coisider  the  case  where  ♦ is  in  the  space  of 

4 by  3 matrices  and  = 1,  = 2,  ^^2  " “ ^13  " ^3  ' ^23  ’ ^23  ' ^2 

V32  ” *•43  “ ^43  “ "I*  *'21  “ ^21  ' *"22  " *^22  “ *'31  * ^31  “ *33  * ^^33 

*41  “ ^41  " *42  ' ^42  ° Hence  the  subset  M of  il>  is  composed  by  the  two 


matrices 


111 


2 2 1 


0 10 


0 0 -1 


and  C 


2 11 


2 2 1 


0 1 0 


0 0-1 


1 


Let  F be  defined  as  the  convex  hull  of  A,  B,  C,  and  D where  A = (0,0,0), 

B “ (3, -3,0),  C = (■'  "I.  0,  ^),  and  D = (5,  - ■^,  1).  Point  D is  in  the  plane 
defined  by  A,  B,  and  C.  It  can  be  shown  that  no  point  in  the  relative 
interior  of  ABCD  is  in  EF(c‘)  or  in  EF(C^),  EF(c')  ACUCD  and  EF(C^)  = AB(^>BD. 
Therefore  EF  = EF(c') nEF(C^)  = {A,d),  l.e.,  EF  has  two  elements,  the  points  A 
and  D.  Moreover,  the  vertices  A and  D are  not  adjacent  in  F. 

Consider  problem  (LMOP) . It  is  well  known  [17]  that  x“eEF(C)  if  and  only 
If  X*  solves  max{XCx:xeF}  for  some  X > 0,  XeR^.  It  is  therefore  natural  to 
consider  systems  (1)  and  (II) , below,  in  order  to  obtain  conditions  for  EF 
to  be  non-empty. 

X^C^  - qI  =0  for  all  1 such  that  C^eM 
n 

X^  > 0 (I) 

X^eR^’,  qeR” 

and 

X^C^  ~ ql  =0  for  all  1 such  that  C^eM 

X^  > 0 (II) 

X^erP,  qeR" 

Clearly,  if  (I)  has  a solution  q®  and  if  the  problem  (Pq®):  max{q®x  : xeF}  has 
an  optimal  solution  then  EF  ^ ({i.  When  F is  bounded  the  condition  on  (Pq®)  is 
always  satisfied  and  can  be  deleted.  Also,  if  for  some  CcM  (or  <]>)  , 

K(C)  ('>{rER”  : Ar=0,  r^O}  # i]),  EF  has  no  elements  since  EF(C)  = ([>.  Whenever 
system  I has  no  solution  it  follows  from  Proposition  A.l  in  [4]  that  the 
closure  of  the  convex  hull  of  K(4>)  contains  a non-trlvial  subspace.  If  in 
addition,  system  II  has  the  unique  solution  q®=0  then,  the  subspace  in  the 
closure  of  the  convex  hull  of  k('!’)  Is  r”.  When  q®  0 is  feasible  in  system  II 
and  the  problem  (Pq®)  ; max{q®x  ; xeF}  has  a unique  optimal  solution  x®,  EF  is 
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non-empty  since  x'eEFCC)  for  all  CeM  . The  two  following  examples  Illustrate 
the  necessity  of  the  condition  of  (Pq°)  to  have  a unique  solution  and  that  EF 
need  not  be  empty  when  the  subspace  in  the  closure  of  the  convex  hull  of 
is  r”.  Let  c’  and  be  the  only  two  elements  of  M where, 


■ 1 

o' 

1 

0 

r-l 

1 

1 

c*  = 

and  C* 

0 

1, 

L 0 ij 

and  let  F be  the  segment  AB  where  A = (-1,0)  and  B = (1,0).  We  have  that 
AeEF(C^),  BeEF(C*),  and  EF(M)  = (|) . Note  that  in  this  case  system  I has  no 
solution,  system  II  has  solution  q“  = (0,1)  and  the  segment  AB  is  optimal 
in  (Pq°).  If  instead  of  the  segment  AB,  F is  defined  as  the  triangle  ABD 
with  D = (0,1)  then,  D becomes  the  unique  solution  of  (Pq“)  and  DeEF.  The 
next  example  Illustrates  the  situation  where  FCR^,  closure  of  the  convex  hull 
of  k(<I>)  Is  but  EF  ^ (}).  Let  C*  = [1,0]  and  = [-1,0]  be  the  only  two 

elements  in  M and  let  F be  the  line  defined  by  the  two  points  A = (0,1)  and 
B = (0,-1).  Thus,  k(C*)  = {psR^  : Pj^>0},  K(C^)  = {pER^  : Pj^<0};  closure  of 
the  convex  hull  of  k(^’)  is  , but  all  points  in  F are  in  EF(C*)  and  EF(C^) 
and  therefore  F = EF. 

3.  Nonnegative  Criteria 

In  this  section  we  concentrate  on  the  case  where  ^ 0 for  all  (l,j). 
Without  loss  of  generality  we  assume  that  for  each  j = {l,2,...,n)  there  is 
a C in  with  the  elements  in  column  n not  all  zeros.  The  next  two  proposi- 
tions suggest  a simple  algoritfira  to  determine  EF^^  when  the  criteria  are  non- 
negative. 

Proposition  3.1:  Let  EF(I^)  be  the  set  of  efficient  points  in  problem 

P(I  ) ; Max{l  X : xeF}.  Then,  EFCEF(I  ). 
n n _ n 
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Propositlon  3.2:  Assume  x“eEF(I  ) and  x®^EF.  Then,  there  is  y“eEF(I  ) such 
* n n 

that  Cy*  ^ Cx"  for  some  C In  t>. 

Proof:  x^^EF  implies  that  there  is  a matrix  Cc<l>  such  that  x°^EF(C).  Let 

e=(l,l, . . . ,l)eR*’  and  consider  the  problem 

(P) : maxfeCx  : Cx^Cx®,  Ax=b,  x^O}. 

If  (P)  is  unbounded  there  is  reR**  satisfying  Cr^O,  Ar=0,  r^O,  eCr>0 

consequently,  Cr^O  and  r^O.  Thus,  r4x“cF  and  r+x®  ^ x®  contradicting  the 

assumption  x®eEF(I  ).  Therefore,  (P)  is  bounded.  Furthermore  y®  solves  (P) 
n 

implies  that  y®eEF(I^)  since  otherwise,  there  would  be  a zeF,  z^®,  Cz^Cy® 
and  eCz>eCy®.  Hence,  any  y®  optimal  in  (P)  is  in  EF(I^)  and  satisfies 
Cy®  > Cx". 

Propositions  3.1  and  3.2  show  that  EF  is  contained  in  the  set  EF(I^) 

and  for  each  x®  efficient  in  and  not  efficient  in  (P)  there  is  an 

element  y®eEF(I  ) and  a matrix  C in  that  can  be  used  to  eliminate  x®  from 
n 

EF(I^) , i.e.,  such  that  Cy®  ^ Cx®.  The  two  propositions  suggest  the  following 
algorithm  to  obtain  the  set  of  efficient  extreme  points  in  (P) . 

Algorithm 

Step  1:  Using  any  algorithm  for  linear  multiple  objective  problems,  determine 

EF  (I  ) . 
ex  n 

Step  2:  For  each  x®eEF  (I  ) and  each  xeEF  (I  ),  x/x® , define  the  matrix 

ex  n ex  n ’ ’ 

C = as  follows: 

C^j  “ U^j  if  Xj-x°j  ^ 0 and 

C.J  - 1..J  it  < 0 fo,  J.1,2 

Note  that  CcM.  If  C(x-x®)  > 0,  x°^EF  and  x®  can  be  discarded  from  EF  (I  ). 

— ex  ex  n 

Otherwise,  choose  another  x from  the  set  of  efficient  points  in  E(Ijj)  and 

repeat  the  procedure.  If  x®  is  not  discarded  after  being  compared  with  the 

elements  of  EF  (I  ) it  follows  from  Proposition  3.2  that  it  is  efficient  in 
ex  n 
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(P) . Proposition  3.1  ensures  that  all  elements  of  can  be  obtained  by 


this  algorithm.  It. is  Interesting  to  note  that  Proposition  3.1  is  valid  when- 
ever the  criteria  are  non-negative.  Therefore  if  for  sensitivity  analysis 


purposes  we  change  the  upper  or  lower  bound  on  any  c^j  we  can  restart  the 


problem  with  the  same  long  as  for  all  (i,j)  and  the  columns 


with  all  components  equal  to  zero  in  all  matrices  in  t remain  unchanged. 
Related  results  were  developed  by  the  author  [3]  for  the  zero-one  linear 
multiple  objective  problem  with  nonnegative  criteria. 


4.  The  General  Linear  Case 


Basically,  the  existing  algorithms  to  determine  the  efficient  extreme 


points  of  (LMOP)  move  from  one  element  of  EF^^(C)  to  an  adjacent  one.  The 


new  extreme  point  is  chosen  by  a test  rule.  The  solution  procedure,  that 


we  propose,  to  obtain  EF  differs  from  the  one  described  in  that  we  solve  an 

ex 


additional  subproblem  to  check  if  the  extreme  point  is  in  To  simplify 


the  notation  we  assume  that  all  extreme  points  of  F,  under  consideration, 
correspond  to  non-degenerate  primal  basis.  At  the  end  of  the  section  it  is 
mentioned  how  this  assumption  can  be  easily  relaxed. 

Initially  we  consider  the  problem  of  how  to  determine  if  a given  extreme 
point  of  F is  or  is  not  efficient  in  (P) . 


B N 

Proposition  4.1;  Let  x*  » (x“  ,x°  ) be  a non-degenerate  extreme  point  of  F. 


Then,  x“eEF^^  if  ad  only  if  the  system 


tN 


(C^  - CgA")p  > 0 


p > 0,  C - (Cg,Cjj)eM 


(4.1) 


has  no  solution. 
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Proof:  Follows  from  the  definitions  of  EF,  M,  and  lemma  2.1  in  [12].  <>■'/  \ 

Proposition  4.1  can  be  written  in  an  operational  form  as:  ”x*eEF^^  if 
and  only  if  the  optimal  value  of  P(x“):  z=max{es:  -1-  I^s  “ 0, 

|i>0,  s^O,  CeM}  is  zero  where,  e = (1,1 l)eR^".  The  best  choice  of 

C •“  (Cp,Cjj)  is  such  that  Cj^  - Ujj  where  U • i«e-,  every  element  of 

should  be  chosen  at  its  upper  limit.  Hence, 

* * * 

Corollary  4.2:  There  is  always  an  optimal  solution  (C  ,p  ,8  ) of  P(x  ) with 

WWA 


The  reader  should  note  that  PCx”)  is  not  a linear  progrannlng  problem 
since  besides  p and  s,  matrix  Cg  is  also  unknown.  Cg  is  p by  m if  A has  full 
row  rank,  and  does  not  depend  on  n.  The  next  propositions  play  an  important 
role  in  the  algorithm  developed  to  solve  P(x®). 

Let  N denote  the  set  of  indices  corresponding  to  the  non-basic  components 
of  x®.  For  each  jeN  let  Cg(j)  be  tbe  p by  m matrix  with  columns  Cg(j),j^ 
defined  as  follows: 


<V.k 

if  A^j  > 0 and 

<V.k 

If  < 0 

k"l,2 , . . . ,m 


where  L « (Lg,Ljj)  and  U “ (Ug,Ujj).  We  refer  to  Cg(j)  as  the  ideal  matrix 
corresponding  to  index  j . It  is  the  p by  m matrix  in  (L-  U„)  that  when  multi- 

D » D 

plied  by  A^j  gives  the  p vector  with  the  lowest  component  values.  Hence, 


Proposition  4.3:  a)  Cg(j).^Ajj  < jeN;  k-1,2.. . . ,m  and  Lg<Cg<Ug, 

b)  Cg(J)  A^j  < CgA^j  for  all  jeN  and  Lg  < Cg  < Ug. 

. * * * * 

Corollary  4.4:  If  P(x  ) has  a feasible  solution  C ■ (Cg,Ujj),  s , p with 

es  >0  then  P(x“,mj^)  defined  as 


V 


max  es 


P(x*,mj) : 


rN 


tN 


....  <0,,.,  . r * x^. 


y > 0,  8 > 0,  aB).^  < (Cg).j  < (Ug).^  1=1,..., 


has  a feasible  solution  with  es  > 0 for  all  0 ^ £ m [for  = 0 (mj^=m)  the 

first  (second)  summation  vanishes  and  P(x®,mj^*m)  « P(x*)]. 

Proof:  By  Proposition  A. 3a)  ^ j ’ **®**“* 


- -4  h-\  <V.i'u-<V.xi“x 


* 
s . 


Therefore  es  > es  >0. 


If  X*  is  not  efficient  in  (P)  it  follows  from  Proposition  2.2  and 
Corollary  A. 2 that  there  is  a matrix  C ■ (C„,C  “U„)eM  and  an  xeF  such  that 
C(x-x*)  > 0.  Thus,  a naive  method  is  to  solve  P(x')  for  all  matrices  Cg 
with  each  column  having  all  elements  at  the  upper  or  lower  bound.  Since 
there  can  be  2”'  of  such  matrices  this  procedure  would  require  a significant 
amount  of  computer  time.  Instead,  we  perform  the  following  implicit  enumera- 
tion which  considers  a sequence  of  linear  programs  in  order  to  solve  the 
non-linear  problem  P(x”). 


Implicit  Enumeration  Algorithm 

Description;  Start  by  solving  P(x®,  mj^=0) . If  v“0  stop,  x®  is  efficient  in  (P) . 
If  v > 0 let  mj“l  and  generate  the  following  two  problems 
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i: 


P(x“,  max  es 


- <V-J  j 


“j  * S* 


and 


y > 0,  s > 0 


P(x®,  max  es 


+ I s 
P 


jeN 

p > 0,  s ^ 0. 

vniere  In  the  notation  PCx®,  mj^=l,t),  t“l  (t=0)  indicates  that  the  column,  in 
C , corresponding  to  m, =1  has  all  its  elements  at  the  upper  (lower)  bound. 

If  the  optimal  value  of  P(x*,  mj^-1,1)  is  zero,  by  Corollary  4.4,  there  is  no 
optimal  matrix  C_  in  P(x®),  with  z > 0 and  having  the  first  column  equal  to 
(U„)  In  this  case  we  do  not  need  to  consider  any  descendant  of  P(x*,  m «1,1) 
and  the  branch  is  fathomed.  If  the  optimal  value  of  P(x*,  m^*!,!)  is  positive 
we  generate  two  new  problems  P(x®,  inj^«2,l,l)  and  P(x*,  mj“2,l,0).  These  two 
problems  are  obtained  by  substituting  Cg(J),2»  P(x*,  ®2^“1»1)»  respectively 

by  (U«)  and  (L„)  ».  We  proceed  in  the  same  way,  i.e.,  branching  on  problems 
with  optimal  value  positive  and  fathoming  those  with  optimal  value  zero  until, 
we  either  conclude  that  x*  is  efficient  in  (P)  or  obtain  a Cj  such  that 
Lg  ^ Cg  ^ Ug  and  z > 0.  An  example  of  a tree  generated  by  the  implicit 
enumeration  algorithm  is  given  in  Figure  1.  In  that  figure,  P(x®,  mj^“3, 1,0,1) 
is  the  problem 

max  es 


U > 0,  s > 0. 
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P(x*,  ini=l,l) 
optimal  value  > 0 


P(x*,  mj^-1,0) 
optimal  value  > 0 


P(x',  mi=2,l,l) 
optimal  value  = 0 
(fathom  branch) 


P(x*,  m -2,1,0) 
optimal  value  > 0 


P(x*,  mj»2,0,l) 
optimal  value  - 0 
(fathom  branch) 


P(x®,  m]^=2,0,0) 
optimal  value  *=  0 
(fathom  branch) 


P(x*,  mj^-3, 1,0,1) 
optiiaal'‘'value  = 0 
(fathom  branch) 


P(x*,  mj-3, 1,0,0) 
optiaal  value  - 0 
(fathom  branch) 


Figure  1;  Example  of  a Tree  Generated  by  the  Implicit  Enumeration  Algorithm. 
In  this  case  x®  is  efficient  in  (P) . 
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The  convergence  of  the  algorithm,  after  solving  a finite  number  of  linear 
programs,  follows  from  Corollary  4.4  and  the  fact  that  the  number  of  matrices 
Cjj  that  can  possibly  be  enumerated  Is  finite. 

In  (LMOP)  the  set  of  efficient  points  is  connected  and  the  algorithms  to 
determine  EF^^(C)  move  from  one  efficient  extreme  point  to  another  by  using 
only  efficient  edges.  In  (P) , the  set  of  efficient  points  Is  not  necessarily 
connected,  as  Illustrated  In  section  2.  However,  whenever  it  is  of  Interest 
to  identify  the  connected  subsets  of  EF  the  following  problem  can  be  solved 
to  determine  If  an  adjacent  edge  to  the  efficient  point  x*  Is  efficient  or  not: 


Q(x®,j):  max  es 

s.t.  (C/-  C^)p  - [CgA^.  - (C^).jlw+  IpS  - 0 

y > 0,  s > 0 


I 


where  Xj  Is  a non-baslc  variable,  in  the  feasible  solution  associated  with  x® 
which  when  introduced  Into  the  basis  generates  an  adjacent  extreme  point  that 
we  denote  by  x^ . The  edge  of  F defined  by  the  pair  (x®,x^)  is  the  segment 
- Ix®,x^]. 


Proposition  4.5:  Let  x®  be  an  efficient  extreme  point  in  F.  Then,  the  edge 
r^  ■ [x*,x^]  is  efficient  In  (P)  if  and  only  If  the  optimal  value  of  Q(x®,j) 

Is  zero. 

Proof:  Follows  by  the  definitions  of  EF,  M,  and  leimna  2.3  in  [12].  ^ 

Corrolary  4.6:  Let  x®  be  an  efficient  extreme  point  in  F.  Then,  the  edge 

r^  > [x*,x^]  is  efficient  in  (P)  if  and  only  if  the  optimal  value  of  the  problem 

problem 


i 


■w 
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Q(x“, j) : max  es 


s.t.  (CgA**  - Uj^)u  - [Cg(j)A”j  - ° 


p>0,  w >0,  w >0,  s^O 


is  zero  where,  Cg(j)  is  the  p by  m matrix  with  columns  defined  as 


S<J) 


•k 


if  A^j  > 0 and 

If  < 0 

Proof:  Follows  from  the  fact  that  the  optimal  values  of  Q(x*,j)  and  Q(x®,j)  are 
equal.  t'/fi'/.  A 

Clearly,  Q(x®,j)  can  also  be  solved  by  the  implicit  enumeration  algorithm. 

To  obtain  we  can  choose  a matrix  C,  for  example  C - y(L+U) , apply 

any  of  the  existing  algorithms  to  solve  (LMOP)  and  at  each  of  its  efficient 

extreme  points  x®  solve  P(x®)  to  determine  if  it  ia  or  not  in  EF^^.  This 

procedure  will  generate  the  whole  set  since  it  is  contained  in  EF^^(C). 

Therefore,  the  method  consists  in  using  the  connected  graph  formed  by  efficient 

extreme  points  and  edges  of  (LMOP)  to  determine  the  corresponding  non-connected 

subgraph  of  (P) . It  is  worthwhile  to  point  out  that  the  algorithms  developed 

by  Yu  and  Zeleny  [36l,  Ecker  and  Kouada  [11],  and  Gal  [16]  to  obtain  efficient 

faces  of  F can  also  be  conveniently  adapted  to  determine  the  efficient  faces  of 

F in  problem  (P) . Similarly  the  Implicit  enumeration  algorithm  can  be  applied 

in  conjunction  with  Steuer's  interval  weights  method  [32].  To  conclude  this 

section,  we  consider  the  non-degeneracy  assumption  made  earlier.  Whenever  x® 

is  a degenerate  basic  feasible  solution  it  is  necessary  to  add  to  problems  P(x®) 

-N 

and  P(x  ,mj)  the  constraints  A*^^p  < 0 for  all  i such  that  the  corresponding 
basic  components  of  x®  have  value  zero. 


j 
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5.  Computational  Results  and  Conclusions 

The  Implicit  enumeration  algorithm  was  tested  by  solving  four  hundred 
problems  of  the  form  of  P(x°).  The  matrices  A were  randomly  generated  In 
the  Interval  [0,20]  with  20%  of  negative  elements  and  the  value  of  the 
right-hand- side  vector  b was  fixed  with  all  components  equal  to  100.  The 
criteria  matrices  were  constructed  as  follows.  For  each  problem  a matrix 
C “ [Cjjl  randomly  generated  in  the  interval  [0,99]  and  the  matrices 
U and  L were  defined  as  U = 6C  and  L = (2-B)C  with  3 * 1.05  or  B = 1.10. 

For  each  problem  an  efficient  point  with  respect  to  C was  generated  by 
solving  the  linear  program  max{XCx  ; Ax=b,  x^O}  with  C as  defined  above  srd 
A » (1,1,1)  or  (1,1,8)  when  p “ 3 and  A = (1,1, 1,1,1)  or  (1,1, 1,1, 6)  when 
p » 5.  The  computer  program  was  written  in  Fortran  and  was  processed  on  the 
mini-computer  Prime  400.  Several  branching  rules  were  tested  in  the 
algorithm  and  are  not  reported  here  due  to  space  limitations.  The  use  of 
different  values  for  A Is  an  attempt  to  Identify  a distinct  behavior  between 
points  that  tend  to  be  In  the  center  of  the  efficient  set,  l.e.  those 
corresonding  to  A*s  with  all  components  equal,  and  those  that  tend  to  be 
In  the  frontier  of  the  efficient  set.  The  computational  results  are  given 
In  Table  5.1.  Each  row  of  the  table  summarizes  the  results  of  a sample  of 
25  problems  P(x“).  m,  n,  p.  A,  6,  y,  and  a denote  respectively  the  number 
of  constraints,  number  of  variables,  number  of  criteria,  vector  of  weights, 
the  number  that  multiply  the  randomly  generated  matrix  C to  determine  !>  “ 
the  mean  c.p.u.  time  to  solve  25  problems  P(x®)  using  the  implicit  enumeration 
algorithm  (this  time  does  not  include  the  determination  of  the  point  x®,  l.e., 
the  c.p.u.  time  to  solve  max{ACx  : Ax”b,  x^O))  and  the  standard  deviation 
of  the  time  to  solve  each  problem  P(x®).  The  last  column  in  Table  5.1  indi- 
cates the  total  number  of  points  that  were  found  efficient  after  solving  the 
25  problems  P(x®).  The  computational  results  show  that  for  B “ 1.05  and 
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Table  5.1 


1 

1 

m ! 

n 

P 

A 

8 

U 

10  ^ sec 

a 

10  ^ sec 

No.  of 
efficient 
points 

10 

20 

3 

(1,1.1) 

1,05 

.17 

.05 

20 

10 

20 

3 

(1.1.8) 

1.05 

.42 

.20 

18 

in 

20 

3 

(1,1,1) 

1.10 

.68 

.29 

15 

in 

20 

3 

(1,1,8) 

1.10 

.32 

.11 

10 

in 

20 

5 

(1,1,1, 1,1) 

1.05 

1.04 

.61 

24 

in 

20 

5 

(1,1,1, 1,6) 

1.05 

2.55 

2.00 

25 

in 

20 

5 

(1,1,1, 1,1) 

1.10 

4.16 

2.01 

20 

in 

20 

5 

(1,1, 1,1, 6) 

1.10 

2.58 

1.19 

15 

in 

50 

3 

(1.1,1) 

1.05 

2.96 

1.81 

18 

in 

50 

3 

(1,1,8) 

1.05 

2.77 

1.21 

18 

10 

50 

3 

(1,1,1) 

1.10 

4.31 

2.35 

15 

in 

50 

3 

(1.1.8) 

1.10 

2.47 

.98 

10 

25 

50 

3 

(1,1,1) 

1.05 

19.90 

14.22 

3 

25 

50 

3 

(1.1,8) 

1.05 

22.67 

14.89 

3 

25 

50 

3 

(1.1.1) 

1.10 

8.01 

5.16 

1 

25 

50 

L — 

3 

(1,1,8) 

1.10 

5.93 

6.70 

0 

m,  n,  and  p fixed  the  mean  computation  time  is  usually  lower  for  equal  weights 
than  for  unequal  weights,  the  opposite  is  true  for  $ = 1.10.  This  fact 
suggests  that  for  small  perturbations  in  the  matrix  C,  i.e.,  for  small  values 
of  6,  problems  PCx*)  for  x“  obtained  with  equal  weights  are  easier  to  solve 
or  equivalently,  the  algorithm  needs  to  explore  fewer  nodes  of  the  tree 
thfln  for  the  case  of  unequal  weights  of  the  type  considered.  It  Is  also 
Interesting  to  note  that  for  vectors  A with  equal  components  the  number  of 
points  found  efficient  tends  to  be  larger  than  for  the  other  weights  consi- 
dered. For  the  ranges  of  n,  m,  and  p tested,  the  mean  computational  time 
seems  to  be  equally  sensitive  to  individual  increases  in  these  parameters. 
However,  the  variance  is  apparently  more  sensitive  to  an  increase  in  the 
number  of  constraints  for  8 “ 1.05  than  for  B “ 1.10. 


Conclusions  and  Topics  for  Further  Research 

We  have  considered  in  this  paper  the  linear  multiple  objective  problem 
with  interval  coefficients,  denoted  by(P).  Properties  of  (P)  were  presented 
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together  with  a problem  to  test  if  a feasible  extreme  point  Is  or  not 
efficient.  This  problem  can  be  used  in  conjunction  with  known  algorithms 
to  determine  the  set  of  efficient  extreme  points  In  (P) . A branch  and  bound 
algorithm  and  computational  results  were  given  and  discussed. 

Potential  extensions  and  topics  for  further  research  are  the  determina- 
tion of  "an  easy  to  check"  necessary  and  sufficient  condition  for  the  exis- 
tence of  an  efficient  point  in  (P)  and  the  consideration  of  situations  where 
the  coefficients  of  the  objective  functions  are  not  independent  random 
variables.  For  the  same  practical  reasons  pointed  out  in  the  introduction 
of  the  paper,  it  is  of  Interest  to  develop  algorithms  to  solve  the  Interval 
coefficient  versions  of  zero-one  linear  multiple  objective  programs,  the 
problem  of  finding  efficient  points  when  the  set  of  discrete  alternatives 
is  given  explicitly  [3]  and  several  other  discrete  programs  suggested  in  [39]. 
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